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In this paper we study the excitation spectrum of graphene in a strong magnetic field, beyond the 
Dirac cone approximation. The dynamical polarizability is obtained using a full 7r-band tight-binding 
model where the effect of the magnetic field is accounted for by means of the Peierls substitution. 
The effect of electron-electron interaction is considered within the random phase approximation, 
from which we obtain the dressed polarization function and the dielectric function. The range of 
validity of the Landau level quantization within the continuum approximation is studied, as well as 
the non-trivial quantization of the spectrum around the Van Hove singularity. We further discuss 
the effect of disorder, which leads to a smearing of the absorption peaks, and temperature, which 
activates additional inter-Landau level transitions induced by the Fermi distribution function. 



I. INTRODUCTION 

One of the most remarkable features of graphene is 
its anomalous quantum Hall effect (QHE), which reveals 
the relativistic character of the low energy carriers in 
this material*^ In fact, the linear electronic dispersion 
of graphene near the neutrality point leads to a rela- 
tivistic quantization of the electrons' kinetic energy into 
non-equidistant Landau levels (LL), with the presence of 
a zero-energy LL, which is the characteristics spectrum 
for systems of massless Dirac fermionsi&i As a conse- 
quence, the excitation spectrum and the screening prop- 
erties in graphene are different from those of a standard 
two-dimensional electron gas (2DEG) with a quadratic 
band dispersion, as it may be seen from the polarization 
and dielectric function in the two cases J±~— 

The Coulomb interaction between electrons in com- 
pletely filled LLs leads to collective excitations and to 
the renormalization of the electronic properties such as 
the band dispersion and the Fermi velocity. These is- 
sues have been studied both theoretically-^— and ex- 
perimentally, in the framework of cyclotron resonance 
experiments ^2r— However, most of the theoretical work 
has been based on the continuum Dirac cone approxima- 
tion, which does not apply when high energy inter-LL 
transitions arc probed^ In recent experimental realiza- 
tion of "artificial graphene" ^ a two-dimensonal nanos- 
tructure that consists of identical potential wells (quan- 
tum dots) arranged in a honeycomb lattice, the lattice 
constant (a ~ 130 nm) is much larger than the one in 
graphene (do ~ 0.142 nm). This provides a way to study 
graphene in the ultra-high magnetic field limit, since a 
perpendicular magnetic field in "artificial graphene" cor- 
rcspondes to an effective field which is (a/ao) 2 ~8x 10 5 
times larger than in graphene. Furthermore, the re- 
cently developed techniques of chemical doping^ and 
electrolytic gating^ have enabled doping graphene with 
ultrahigh carrier densities, where the band structure is 
no longer Dirac-like and one should take into account 
the full 7r-band structure including the Van Hove singu- 



larities (VHS). 

In this paper, we present a complete theoretical study 
of the density of states (DOS), the polarizability and di- 
electric function of graphene in a strong magnetic field, 
calculated from a 7r-band tight-binding model. The mag- 
netic field has been introduced by means of a Peierls 
phasej&2& and the effect of long range Coulomb interac- 
tion is accounted for within the random phase approxi- 
mation (RPA). Our method allows us to study the effect 
of temperature, which leads to the activation of addi- 
tional inter-LL transitions. We also study the effect of 
disorder in the spectrum, which leads to a smearing of 
the resonance peaks. 



II. DESCRIPTION OF THE METHOD 

In this section we summarize the method used in the 
numerical calculation of the polarizability of graphene in 
the QHE regime4£ A monolayer of graphene consists of 
two triangular sublattices of carbon atoms with an inter- 
atomic distance of a w 1.42 A. By considering only first 
neighbor hopping between the p z orbitals, the 7r-band 
tight-binding Hamiltonian of a graphene layer is given 
by 
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where at (hi) creates (annihilates) an electron on sublat- 
tice A (B) of the graphene layer, and ttj is the nearest 
neighbor hopping parameter, which oscillates around its 
mean value t«3 eV<2I The second term of H accounts 
for the effect of an on-site potential Vi, where n.j = c\ci is 
the occupation number operator. For simplicity, we omit 
the spin degree of freedom in Eq. ([1}, which contributes 
only through a degeneracy factor. In our numerical calcu- 
lations, we use periodic boundary conditions. The effect 
of a perpendicular magnetic field B = Bz is accounted 
by means of the Peierls substitution, which transforms 



2 



the hopping parameters according 
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where $0 = hc/e is the flux quantum and A is the vector 
potential, e. g. in the Landau gauge A = (—By ,0,0). 
We will calculate the DOS and the polarization function 
of the system by using an algorithm based on the evo- 
lution of the time-dependent Schrodinger equation. For 
this we will use a random superposition of all basis states 
as an initial state \ip) (see e. g. Refs|28||2£" 
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where are random complex numbers normalized as 
Si \ a i\ = 1- The DOS, which describes the number 
of states at a given energy level, is then calculated as a 
Fourier transform of the time-dependent correlation func- 
tions 
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with the same initial state defined in Eq. ([3]). The dy- 
namical polarization function can be obtained from the 
Kubo formula 30 



n(q,w) = - J rfre" 



([p(q 1 r),p(-q > 0)]>, (5) 



where V denotes the volume (or area in 2D) of the unit 
cell, p (q) is the density operator given by 



P(q) =J2c\c l exp (iq-rj, 



(6) 



and the average is taken over the canonical ensemble. For 
the case of the single-particle Hamiltonian, Eq. (J5j) can 
be written as2^ 
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is the Fermi-Dirac distribution operator, f3 = 1/UbT 
where T is the temperature and ks is the Boltzmann con- 
stant, and p is the chemical potential. In the numerical 
simulations, we use units such that H = 1, and the aver- 
age in Eq. (|7|) is performed over the random superposition 
Eq. ((3J). The Fermi-Dirac distribution operator n F (H) 
and the time evolution operator e~ tHT can be obtained 
by the standard Chcbyshev polynomial decomposition^ 




FIG. 1: Constant energy contours obtained from the band 
dispersion Eq. ()13[) . The thick black lines correspond to dis- 
persion at the VHS |e| = t. Notice that the CEC are cen- 
tered around the Dirac points for |e| < t and around T for 
|e| > t. For illustrative reasons, the hexagonal BZ is shown 
in white. For undoped graphene, the valence and conduction 
bands touch each other at the vertices of the hexagon, the so 
called Dirac points (K and K'). The Van Hove singularity lies 
at the M point, and we have defined 8 as the angle between 
the wave-vector q and the fc^-axis. 



Long range Coulomb interaction is considered in the 
RPA, leading to a dressed particle-hole polarization 
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is the Fourier component of the Coulomb interaction in 
two dimensions, in terms of the background dielectric 
constant k. Furthermore, the dielectric function of the 
system is calculated as 



e(q,w) = l-V(g)II(q,w). 



(11) 



The collective modes lead to zeroes of e(q, w), and their 
dispersion relation is defined from 



Re e(q,uj p i) = I - V (q)TL(q, uj pl ) = 0, 



(12) 



which leads to poles in the response function The 
technicalities about the accuracy of the numerical results 
have been discussed elsewherei 28 i 29 ' 31 Here we just men- 
tion that the efficiency of the method is mainly deter- 
mined by three factors: the time interval of the propa- 
gation, the total number of time steps, and the size of 
the sample. The method is more efficient in the presence 
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FIG. 2: Left: DOS of a monolayer of graphene in a magnetic field, for two different values of B. Right: zoom of the low energy 
region of the spectrum. The red vertical lines indicate the position of the LLs in the continuum Dirac cone approximation, Eq. 
(|I9|1 , We have used a sample made of 4096x4096 atoms. 



of strong magnetic fields. Because for weak fields (e.g., 
B < 1 T) the energy difference between LLs become very 
small, this makes that the total number of time steps and 
the size of the sample have to be large enough to provide 
the necessary energy resolution in the numerical simula- 
tion. 



III. DENSITY OF STATES AND EXCITATION 
SPECTRUM 

In this section we study the DOS and the excitation 
spectrum of a graphene layer in a magnetic field, ne- 
glecting the effect of disorder and electron-electron inter- 
action. The B = dispersion relation of the tt bands ob- 
tained from a tight-binding model with nearest-neighbor 
hopping between the p z orbitals is 



e(k) = At^kl 
where A = ±1 is the band index and 

k = l + 2e^ a / 2 cos( ^ 
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The band dispersion Eq. (fl~3"|) consists of two bands that 
touch each other in the vertices of the hexagonal Brilloin 
zone (BZ) (Fig. [TJ, which are the so called Dirac points. 
In the absence of longer range hopping terms, the band 



structure is electron-hole symmetric, and the constant 
energy contours (CEC) obtained from Eq. (| 13[) are shown 
in Fig. [T] For undoped graphene (/x = 0), which is the 
band filling that we will consider all along this paper, the 
Fermi surface consists of just six points at the vertices 
of the BZ. In this case, the low energy excitations can 
be described by means of an effective theory obtained 
from an expansion of the dispersion Eq. (|13[) around 
the K points. This leads to an approximate dispersion 
e(k) ~ Aupfc, where vf = Zta/2 is the Fermi velocity. 

If we now consider the effect of a perpendicular mag- 
netic field, the Landau quantization of the kinetic energy 
leads to a set of LLs, which can be described from the 
semiclassical conditional^ 
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is the area enclosed by the cyclotron orbit C in mo- 
mentum space [for circular orbits S(C) is just 7rfc 2 ], 
I b = \ffrc/eB is the magnetic length, n is the LL index 
and r(C) is the Berry phase. In graphene, as we will dis- 
cuss in the next section, T(C) = tt for orbits around the 
K and K' points, and for orbits around the T point^ 
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FIG. 3: DOS obtained from Eq. QJ at B = 500T. The insets 
show a zoom of the DOS around the VHS |e| ~ t and at 
|e| w 2.5t. 

From the energy dependence of S(C) one can calculate 
the energy of the Landau levels as 

e - = * r + 2-^rJJ' (17) 

where S f_1 (x) is the inverse function to S(x). Using Eq. 
(|T7|) . it is easy to check that the LL quantization corre- 
sponding to a low energy parabolic band e(k) = k 2 /2mt, 
(where mj, is the effective mass) with a r(C) = Berry 
phase, is e„ = w c (n + 1/2), where uj c = eB/nif, is the 
cyclotron frequency. On the other hand, a linearly dis- 
persing band as the one for graphene leads to a LL quan- 
tization around the Dirac points as 

i — i 

e\,n = Ae„ = A- — v2n oc V Bn. (18) 



A. Density of states 

The DOS close to the Dirac point can be approximated 
bySi 

d Dira c(e)«— 4 (19) 

where A c = 3v / 3a 2 /2 is the unit cell area. In Fig. [2] 
we show the DOS for two different values of the mag- 
netic field. The black line corresponds to the numerical 
tight-binding result obtained from Eq. Near e = 
we notice the presence of a zero energy LL surrounded 
by a set of LLs whose separation decreases as the en- 
ergy increases, leading to a stacking of the LLs as we 
move away from the Dirac points. The presence of a 
finite broadening in these LLs is due to the energy reso- 
lution of the numerical simulations, which is limited by 



the number of atoms used in the calculation, as well as 
the total number of time steps, which determines the ac- 
curacy of the energy eigenvalues. In order to check the 
range of validity of the continuum approximation, in the 
right-hand side of Fig. [5] we show a zoom of the pos- 
itive low energy part of the DOS for the two values of 
B, comparing the DOS obtained with the full 7r-band 
tight-binding model Eq. (0} [black lines] to the Dirac 
cone approximation of Eq. (|19p [vertical red lines] . Con- 
trary to multi-layer graphenes, for which trigonal warp- 
ing effects are important at rather low energies^ we see 
that the deviations of the LL positions in the contin- 
uum approximation Eq. (TPI)]) with respect to the full 
7r-band model are weak even at energies of the order 
of e ~ 0.3t ~ leV, in agreement with magneto-optical 
transmission spectroscopy experiments^ 

A much less investigated issue is the effect of the mag- 
netic field on the DOS around the VHS |e| « t. For illus- 
trative reasons we show in Fig. [3]the numerical results for 
the DOS of a graphene layer at an extremely high mag- 
netic field of B = 500T42 At this energy the LL quan- 
tization is highly nontrivial because of the saddle point 
in the band structure at which there is a transition from 
CECs encircling the Dirac points, to CECs encircling the 
r point, as it can be seen in Fig. [TJ Because in the semi- 
classical limit, the cyclotron orbits in reciprocal space 
follow the CECs, we have that at the saddle point there 
is a change in the topological Berry phase T(C) from 
r(C) = ±7t for orbits encircling the Dirac points (|e| < t) 
to r(C) = for orbits encircling the T point (|e| > t)M- 
The different character of the cyclotron orbits at both 
sides of the saddle point leads to two series of LLs, with 
different cyclotron frequencies uj c = eB/nib, that merge 
at the VHS, as it may be seen in the left-hand side inset 
of Fig. [3J Because of the effective mass mj below the 
VHS is larger than the one above it (the band below the 
saddle point is flatter than above it), the cyclotron fre- 
quencies are also different w c (|e| < t) < uj c (\e\ > t), and 
consequently the LLs are more separated for |e| > t than 
for |e| < t. The possibility of placing the chemical poten- 
tial at the VHS would bring the chance of studying highly 
anomalous inter-LL transitions, due to the different sep- 
aration of the LLs above and below the VHS. Finally, at 
an even higher energy, the LL quantization is quite sim- 
ilar to that of a 2DEG with a parabolic dispersion, with 
a set of roughly equidistant LLs, as it may be seen in the 
right-hand side inset of Fig. [3] for |e| rj 2.5t. However, 
we emphasize that for realistic values of magnetic field, 
the LL quantization in graphene is inappreciable in this 
range of energies, and the DOS for energies |e| > 0.7t is 
similar to the DOS at B = 0,™ as it may be seen in Fig. 

m 



B. Particle-hole excitation spectrum 

The particle-hole excitation spectrum (PHES) for non- 
interacting electrons, which is the part of the lu — q 
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FIG. 4: — ImII(q, uj) for different values of wave-vector q and strength of the magnetic field B. The angle 8 defines the 
orientation of the wave- vector in the Brilloin zone (see Fig. [TJ. Plots (a) and (c) show the polarization in the whole energy 
range. For comparison we show the polarization at B — 0. In plots (b) and (c) we show the low energy part of the spectrum. 
The vertical black lines signal the energy of the particle-hole processes expected from Eq. (|20|l . 



plane where Imll(q, uS) is non-zero, defines the region of 
the energy-momentum space where particle-hole excita- 
tions are allowed. For undoped graphene (fi = 0), the 
particle-hole excitations correspond to inter-band tran- 
sitions across the Dirac points. In Fig. 2] we show 
— ImII(q, uj) for different values of wave- vector and mag- 
netic field. Two different orientations of q are shown, 
namely along the T-M and T-K directions. First, due to 
the low energy linear dispersion relation and to the ef- 
fect of the chirality factor (or wave-function overlap) that 
suppresses backscattering, we observe that the strongest 
contribution to the polarization is concentrated around 
ui w vpq. In fact, at B = it was shown^ that 
ImII(q, a;) ~ q 2 (uj 2 — v 2 F q 2 )~ x l 2 , which implies an infi- 
nite response of relativistic non-interacting electrons in 
graphene at the threshold uj = v_f<?- The main difference 
of the PHES at finite magnetic field with respect to its 
-6 = counterpart in this low energy range, is that in the 
B case II(q, uj) presents a series of peaks of strong 
spectral weight, due to the LL quantization of the kinetic 
energy, that we will discuss in more detail below. We no- 
tice that, for the realistic values of magnetic field used in 
Fig. HI the spectrum at finite magnetic field roughly co- 



incides with the the one at B = at high energies. This 
is due to the almost negligible effect of the magnetic field 
on the DOS at energies |e| > 0.7t for B < 50T, as we 
saw in the previous section. This part of the spectrum 
is dominated, as in the B = case ) 31 i 36 by a peak of 
ImII(q, uj) around ui ~ 2t, which is due to particle-hole 
transitions between states of the Van Hove singularities 
of the valence and the conduction bands at e rj — t and 
e ~ t respectively. 

However, the low energy part of the spectrum is com- 
pletely different to its zero magnetic field counterpart, 
and it is dominated by a series of resonance peaks at 
some given energy values. For the undoped case stud- 
ied here, the possible excitations correspond to inter-LL 
transitions of energy uj n _ n i ~ e n i + e„, where n! is the LL 
index of the particle in the conduction band (A = +1) 
and n is the LL index of the hole in the valence band 
(A = —1). In the continuum approximation, they have 
an energy 

= V2{v F /l B )(Vri + y/n). (20) 

The energy corresponding to each of these transitions 
is indicated by a black vertical line in Fig. HJb) and 
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(d), where we show a zoom of — ImII(q, uj) that corre- 
sponds to the low energy LL transitions about the Dirac 
points. Notice that, in contrast to a standard 2DEG with 
a parabolic dispersion and equidistant LLs, the relativis- 
tic quantization of the energy band in graphene makes 
that in a fixed energy window at high energies, there are 
more possible inter-LL excitations from the level n in the 
valence band to the level n' in the conduction band, than 
at low energies. As a consequence, there is a stacking of 
neighboring LL transitions as we increase the energy of 
the excitations, which manifests itself in a continuum of 
possible inter-LL transitions from a given energy, which 
for the case of Fig. gfb) at B = 20T is u>/t > 0.25. 

Contrary to what one can naively expect, only at very 
strong magnetic fields and for small values of w and q, 
the peaks of — ImII(q, cj) occur at the energies given by 
Eq. ([2"U]) . In fact, we can see in Fig. EJd) that for 
B = 50T and for the smaller value of q shown (red line) , 
the peaks of — ImII(q, w) match very well the energies for 
the inter-LL transitions given by Eq. (|20[) . However, at 
weaker magnetic fields and/or larger wave- vectors, the 
peaks of the polarization function do not coincide any 
more with every of the inter-LL transitions given by Eq. 
([2Uf . In fact, there is a series of peaks of ImII(q, uj) , corre- 
sponding to regions of the PHES of high spectral weight, 
which can be understood from the form of the wavefunc- 
tions of the electron and the hole that overlap to form 
an electron-hole pair. A detailed discussion about the 
structure of the PHES in graphene in comparison with a 
2DEG can be found in Ref. |8j. Here we just remember 
that the modulus of the LL wavefunction, due to the ze- 
ros of Laguerre polynomials, presents a number of nodes 
that depend on the LL index n. On the other hand, 
the existence of an electron-hole pair will be possible if 
there is a finite overlap of the electron and hole wavcfunc- 
tions, which will define the form factor for graphene in 
the QHE regime, J- n .n'{^- Because the node structure 
of the single-particle wavefunctions will be transfered to 
\Fn,n' (<l)| 2 1 an together will lead to a highly modulated 
spectral weight in the PHES, as it is seen in Fig. E£b) 
and (d). 



IV. COLLECTIVE MODES 

In the previous section we have discussed the excita- 
tion spectrum in the absence of electron-electron inter- 
action. In this section we include in the problem the 
effect of long range Coulomb interaction. The polariza- 
tion and dielectric functions are calculated within the 
RPA, Eqs. © and (fTTJ). Within this framework, the ex- 
istence of collective excitations will be identified by the 
zeros of the dielectric function or equivalcntly by the di- 
vergences of the loss function — Im[l/e(q, w)], which is 
proportional to the spectrum measured by Electron En- 
ergy Loss Spectroscopy (EELS)i2I In Fig. [5] we show the 
loss function of graphene in a magnetic field, as compared 
to the one at B = 0. At B = the main structure is the 
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FIG. 5: (a) Loss function — Im l/e(q, w) in the RPA, for 
different values of wave-vector. The results for graphene in a 
magnetic field of B — 50T are compared to the B — case, 
(b) Zoom of the low energy part of the spectrum. 



broad peak at w ~ 2t, associated to the 7r-plasmon^i 
When the graphene layer is subjected to a strong per- 
pendicular magnetic field, Im[l/£(q,w)] presents a series 
of prominent peaks at low frequencies, associated to col- 
lective modes in the QHE regime, as it may be seen in 
Fig. \5\b). Similarly to the single-particle case discussed 
in Sec. IIII B[ the strength of the peak is determined by 
the Coulomb matrix elements V(q)\J- n , n >(q)\ 2 , which de- 
pends strongly on the wave- vector q. We emphasize that 
these modes can not be understood as a simple many- 
body renormalization of the dispersionless inter-LL tran- 
sitions given in Eq. (|2"tJj) . because only the low energy 
and long wavelength modes have their non-interacting 
counterpart ui n . n > associated to a specific single-particle 
electron-hole transition with well defined indices n and 
n'. As we go to higher energies and/or weaker mag- 
netic fields, the relativistic LL quantization of graphene 
leads to a so strong LL mixing that the collective modes 
cannot be labeled any more in terms of single-particle 
excitations ; 38 ' 39 as in the case of a 2DEG with a quadratic 
dispersion and a set of equidistant LLs4S 

In Fig. [6] we compare the real part of the dielectric 
function Re e(q, ui) for zero and finite magnetic field. 
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FIG. 6: Re e(q,cj) for B = (red lines) and B = 50T (black 
lines), for the values of the wave- vectors used in Fig. 



The zeros of Re e(q, u>) correspond to the frequencies of 
the undamped collective modes. Contrary to the B = 
case, for which there is no collective modes for undoped 
graphene at the RPA level, at B ^ we observe a num- 
ber of well defined zeros for Re e, which correspond to 
coherent and long-lived linear magneto-plasmonsJ^ The 
divergence of ImII(q, lj) at uj = vpq and the absence of 
backscattering in graphene make that, as we increase the 
wave-vector q, the more coherent collective modes are 
defined also for higher frequencies, namely around the 
threshold u ps Vpq, which is the frequency associated to 
the highest peaks in Fig. [5jb) . For even higher energies, 
the main contribution to the modes of large frequencies 
are inter-LL transitions between well separated LLs, with 
the subsequent reduction in the overlap in the electron- 
hole wavefunctions. Therefore, the collective modes will 
suffer a stronger Landau damping as higher frequencies 
are probed. 



A. Effect of disorder 

Now we focus our attention on the effect of disorder on 
the DOS and on the excitation spectrum of graphene in 
the QHE regime. In general, disorder leads to a broad- 
ening of the LLs, with extended (delocalized) states near 
the center of the original LL, and localized states in the 
tails. We consider here two different kinds of disorder, 
namely random local change of the on-site potentials Vi, 
which acts as a chemical potential shift for the Dirac 
fcrmions, and random renormalization of the hopping 
amplitudes ty, due e. g. to changes of distances or an- 
gles between the carbon p z orbitals. They enter in the 
single-particle Hamiltonian as given in Eq. The ef- 
fect of correlated long-range hopping disorder has been 
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FIG. 7: (a) DOS of clean graphene (black lines) and of dis- 
ordered graphene with a random on-site potential (red lines) 
and with a random renormalization of the hopping integrals 
(blue lines). The inset shows the smearing of the VHS peak 
due to disorder, (b) Zoom of the low energy part of the spec- 
trum. 



shown to lead to a splitting of the n — LL , ' origi- 
nated from the breaking of the sublattice and valley de- 
generacy. Other kinds of disorder as vacancies create 
midgap states that make the n = LL to remain robust, 
whereas the rest of LLs are smeared out due to the effect 
of disorder^ In Fig. [7] we show the DOS of graphene in 
a perpendicular magnetic field of B — 50T for different 
kind of disorder, as compared to the clean case. We let 
the on-site potential «j to be randomly distributed (in- 
dependently on each site i) between —v r and +v r . On 
the other hand, the nearest-neighbor hopping is ran- 
dom and uniformly distributed (independently on sites 
between t — t r and t + t r . At high energies, as we 
have seen in Sec. IIII A[ the DOS for this strength of 
the magnetic field is quite similar to the DOS at B = 0. 
Therefore, as in the zero field case ; 29 i 31 the main effect 
away the Dirac point is a smearing of the VHS at |e| = t, 
as it is observed in the inset of Fig. [7f a). 

Well defined LLs occur around the Dirac point, and the 
effect of disorder on the peaks is observed in Fig. [JJ where 
we show a zoom of the low energy part of the DOS around 
e = 0. Both kinds of disorder (random on-site potential 
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FIG. 8: (a)-(c) Non-interacting single-particle excitation spectrum of graphene in a magnetic field of B = 50T. as defined by 
— ImII(q, u), for different kinds of disorder and for different values of the wave- vector q. The black vertical lines signal the 
energy of the electron-hole processes defined by Eq. ((20]). (d)-(f) Loss function — Im[l/e(q, uj)] in the RPA. (g)-(i) Real part of 
the dielectric function Re e(q, cj) in the RPA. 



and random hopping) leads to a similar effect, producing 
a broadening of the LLs. We can also observe, especially 
for the highest LLs shown in Fig. E[b), a tiny but ap- 
preciable redshift of the position of the center of the LLs 
with respect to its original position, in agreement with 
previous worksi 42 ' 43 The full self-consistent Born approx- 
imation calculations for graphene with unitary scatterers 
of Ref. HH leaded to a rather significant change in the po- 
sition of the LLs towards higher energies. However, the 
exact transfer matrix and diagonalization calculations of 
Ref. |45| found only a small shift of the LL position for 
very strong disorder. 

We now study the effect of disorder on the PHES. In 
Fig. H{a)-(c) we show — ImII(q, ui) for different values 
of q, and for different kinds of disorder. First we no- 
tice a smearing of the resonance peaks associated to the 
LL broadening due to disorder. Whereas for low ener- 



gies the position of the resonance peaks of disordered 
graphene coincides with the position for the clean case, 
we observe a redshift of the resonance peaks as we con- 
sider inter-LL transitions of higher energies. This is due 
to the change in the position of the high energy LLs of 
disordered graphene with respect to the original LLs, as 
we have discussed previously [see Fig. Efb)]. 

The presence of disorder will also affect the dispersion 
relation and coherence of the collective modes when the 
effect of electron-electron interaction is included. In par- 
ticular, the effect of a short range disorder on graphene in 
a magnetic field can lead, due to the possibility of inter- 
valley processes associated to the breakdown of sublattice 
symmetry, to the localization of some collective modes on 
the impurity 4& In Fig. [5fd)-(f) we show the effect of a 
random on-site potential or a random hopping renormal- 
ization on the loss function. First, we can observe an 
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FIG. 9: DOS of clean graphene for different values of temper- 
ature T and magnetic field B. The shaded area is a sketch of 
the Fermi-Dirac distribution function for each case. 



important attenuation of the intensity peaks due to dis- 
order. Second, as we have discussed above, the position 
of the peaks of disordered and clean graphene coincides 
at low energies, but not at high energies, where the res- 
onance peaks of the loss function of disordered graphene 
are shifted with respect to clean graphene. Although we 
obtain a similar renormalization of the spectrum for the 
two kinds of disorder considered here, we reiterate that 
this effect is highly dependent on the type disorder con- 
sidered, as well as the theoretical method used to obtain 
the spectrum. 

Finally, we mention that the disorder LL broaden- 
ing leads to an amplification of the LL mixing discussed 
above. As a consequence, some collective modes which 
are undamped for clean graphene, start to be Landau 
damped due to the effect of disorder. Indeed, in Fig. 
[HKg)- (i) we can see how Eq. (fT2)) , which is the condi- 
tion for the existence of coherent collective modes, is 
fulfilled more times for the clean case than for the dis- 
ordered membranes, for which the collective modes are 
more highly damped. 



B. Effect of temperature: thermally activated 
electron-hole transitions 



In this subsection we discuss the effect of temperature 
on the excitation spectrum, which enters in our calcu- 
lations trough the Fermi-Dirac distribution function Eq. 
([8]) . Temperature can activate additional inter-LL transi- 
tions due to free thermally induced electrons and holes in 
the sample. In Fig. |H]we sketch the LLs of undoped clean 
graphene that are thermally activated, for two different 
strengths of the magnetic field. For this, the correspond- 
ing Fermi-Dirac distribution function is sketched on each 
plot, indicating the LLs that can be partially populated 



with electrons (holes) in the conduction (valence) band. 
Notice that at T = and for fi = 0, np(e) is just a step 
function that traverses the n = LL. Of course, the num- 
ber of activated LLs, which are those crossed by the tail 
of riF(e), grows as we increase the temperature and/or as 
we decrease the magnetic field. The population effect due 
to the thermal excitations of carriers has been observed 
by far infrared transmission experiments^! 

In Fig. [10] we show the single-particle polarization 
and the loss function for two values of temperature and 
magnetic field. At room temperature and for the rather 
strong magnetic fields considered in our calculation, the 
allowed electron-hole transitions are the same as in the 
zero temperature limit (see the top panels of Fig. [5]). 
Therefore, the peaks of ImlT for T = 300 K are cen- 
tered at the frequencies of inter-LL transitions marked 
by the black vertical lines, which accounts only for the 
usual inter-band transitions across the Dirac point. For a 
considerable higher temperature of T = 10 3 K there are 
additional electron-hole transitions (some of them intra- 
band processes, especially important at low frequencies) 
which are now allowed due to the effect of temperature, as 
marked by the red vertical lines in Fig. fTUT aVfb). These 
thermally activated inter-LL transitions at high temper- 
atures contribute to the additional spectral weight of the 
PHES of Fig. Q3n»-(b). Finally, in Fig. HDtc)-(d) we 
show the loss functions corresponding to the magnetic 
fields and temperatures discussed above. As we have dis- 
cussed above, the peaks of Im 1/e correspond to the posi- 
tion of collective excitations. Here we find, in agreement 
with previous tight-binding and band-like matrix numer- 
ical methods^ a weak but appreciable renormalization 
of the collective mode peak position as a function of tem- 
perature. This temperature dependence of the collective 
mode is easily noticed by comparing the red and blue 
peaks at u> « 0.12* of Fig. HOTd). 



V. CONCLUSIONS 

In conclusion, we have studied the excitation spectrum 
of a graphene layer in the presence of a strong magnetic 
field, using a full 7r-band tight-binding model. The mag- 
netic field has been introduced by means of a Peierls sub- 
stitution, and the effect of long range Coulomb interac- 
tion has been considered within the RPA. For realistic 
values of the magnetic field, the LL quantization leads 
to well denned LLs around the Dirac point, whereas the 
DOS at higher energies is rater similar to the one at zero 
field. However, we have shown that in the ultra-high 
magnetic field limit^ for which the magnetic length is 
comparable to the lattice spacing, the LL quantization 
around the Van Hove singularity is highly nontrivial, with 
two different sets of LLs that merge at the saddle point. 

Our results for the polarization function shows that, at 
high energies, the PHES is dominated, as in the B = 
case,— by the 7r-plasmon, which is associated to the en- 
hanced DOS at the VHS of the tt bands. The low energy 
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FIG. 10: (a)-(b) — ImII(q, oj) for two different values of T and B. The black vertical lines denote the position of the inter-LL 
transitions given by Eq. (|20|l which are possible at T = 0. The red vertical lines correspond to the energy of thermally activated 
inter-LL transitions, (c)-(d) Loss function — Im[l/e(q, u>)] for the same values of T and B. 



part of the spectrum is however completely different to 
its zero field counterpart. The rclativistic LL quantiza- 
tion of the spectrum into non-equidistant LLs lead to a 
peculiar excitation spectrum with a strong modulation 
of the spectral weight, which can be understood in terms 
of the node structure of the electron-hole wavefunction 
overlap £ Furthermore, we have shown that the presence 
of disorder in the sample lead to a smearing of the reso- 
nance peaks of the loss function, and to an enhancement 
of the Landau damping of the collective modes. Finally, 
we have studied the effect of temperature on the spec- 
trum, and shown that it can activate additional inter-LL 



transitions, effect which is especially relevant at low mag- 
netic fields. 
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